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Abstract 

Ph! A decay of weakly metastable phase coupled to two-dimensional Liouville gravity is con- 

sidered in the semiclassical approximation. The process is governed by the "critical swelling" , 
where the droplet fluctuation favors a gravitational inflation inside the region of lower energy 



phase. This geometrical effect modifies the standard exponential suppression of the decay 
rate, substituting it with a power one, with the exponent becoming very large in the semi- 
classical regime. This result is compared with the power-like behavior of the discontinuity 
in the specific energy of the dynamical lattice Ising model. The last problem is far from 
being semiclassical, and the corresponding exponent was found to be 3/2. This exponent is 
expected to govern any gravitational decay into a vacuum without massless excitations. We 
conjecture also an exact relation between the exponent in this power-law suppression and 
the central charge of the stable phase. 

Preliminaries. Nucleation mechanism in the first order phase transitions is relatively 
well understood. Nucleation is the main channel of a metastable phase decay in the case of 
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weak metastability (see e.g., [1]). The decay rate turns out to be exponentially small in the 
energy gap between the metastable and stable phases. The leading exponential estimate is 
essentially simple, controlled by the energy of the "critical droplet" . It will be schematically 
repeated (in the particular case of two-dimensional space) in the next section. More sophis- 
ticated calculation allows for evaluation of the pre-exponential factor as well [2]. Somewhat 
later similar mechanism has been studied in [3] (and then independently in [4]) in the con- 
text of quantum field theory, the thermal fluctuations responsible for the formation of critical 
droplet being substituted now by the quantum tunneling. The process is commonly treated 
dynamically in this case, the role of the critical droplet being played by the "materialization 
of a bubble of true vacuum" [4]. In the euclidean version the critical droplet is a kind of 
"instanton", in the sense that it is a localized solution of the classical equations of motion 
(being however an unstable extremum of the action). 

Another facet of this phenomenon can 



"false" vacuum 




be seen when looking at the specific free 
energy as an analytic function of the 
energy gap /x between the phases (e.g., 
/j is proportional to the magnetic field 
strength h in the case of low temperature 
phase ferromagnet in a weak field) . When 
continued to the values corresponding to 
the metastable state, the free energy ex- 
hibits a branch cut with pure imaginary 
discontinuity, directly traceable to the 
"critical droplet" fluctuations (in partic- 
ular, the discontinuity is exponentially 
small in the absolute value of fj,, as in the 
Eq.© below) [5-7]. It is essential that 
the singularity appears only in the ther- 
modynamic limit, when the size of the 
system is sent to infinity. In the field- 
theoretic context the above specific free 
energy is interpreted as the vacuum en- 
ergy density, and the presence of the imaginary part simply means that the "false vacuum" 
is the quantum-mechanically metastable state, thus suggesting direct interpretation of the 
above discontinuity in terms of the decay probability [3,4]. In the context of the thermal 
phase transition, a close qualitative correlation between the discontinuity of the free energy 
and the rate of decay of the metastable phase is widely assumed, while exact relation has 
been established in the limit of weak metastability only [2]. 

The subject attracts much attention. Whereas the stat-mechanical aspects are of obvious 
interest in a wide range of disciplines from chemistry and metallurgy to astrophysics, the 
interest in the "false vacuum" decay is mainly due to the natural concern about the fate of 



Figure 1: A critical droplet fluctuation inside the 
spatially infinite metastable phase 
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our Universe, if it faces the vacuum metastability disaster [8] 3 . There are even more worries 
after Coleman and De Luccia had studied the effects of gravity on this decay [9] . Although 
the Einstein gravity hardly influences the decay rate, the "realistic" parameters taken into 
account, the vacuum formed "a microsecond after" hardly can comfort human beings. For 
more similar horrors and related ethical problems see [10], and if even that is not enough, 
also [11]. Another interesting, although less dramatic, are the implications in the cosmology, 
the scenario where the whole drama is well in the past [12]. 

Turn to less apocalyptic matters. It has been observed recently that the effect of quantum 
gravity in two dimensions changes qualitatively the character of the discontinuity in the 
specific free energy [13]. A lattice version of quantum gravity coupled to Ising spins, known 
as the dynamical lattice Ising model, is exactly solvable [14] through the matrix model 
technique. The vacuum energy as a function of the temperature and magnetic field allows 
analytic expression. In the case of ordered phase, at small absolute values of the magnetic 
field there is indeed a branch cut, the discontinuity being suppressed as a power rather then 
exponentially. 

Some details of this observation are briefly 
recapitulated in one of the subsequent sec- 
tions. It suggests naturally that this power 
law also governs the decay rate in the limit 
of weak metastability. We find it therefore 
instructive to repeat the semiclassical critical 
droplet calculation taking into account rel- 
evant effects of the quantum Liouville grav- 
ity. The latter, rather then the Einstein grav- 
ity, governs the dynamics of the geometry in 
two dimensions [15]. The Liouville gravity 
is characterized by a dimensionless coupling 
constant, related to the central charge c of 
the massless modes around the vacuum. The 
semiclassical limit corresponds to large nega- 
tive c — > — oo. In this limit we find that the 
decay is controlled by a gravitational version 
of critical droplet, which we call the "criti- 
cal swelling" . This localized distortion of the 
metric is induced by the nucleus of the "true" 
vacuum. It looks more like a "cap" at the ex- 
treme semiclassical limit c/i — > — oo and like 
a "bubble" in the opposite weak metastabil- 
ity limit cfi — > —0. In the last case the decay 
rate behaves as /i -c / 6 at fi — > 0, the power- 
law replacing the exponential suppression in the flat. 




(a) 



'bubble" 




(b) 



Figure 2: "Critical swellings in the case s < 1 
(a) and s > 1 (b). 
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This result is semiclassical. It is natural to expect in general some power behavior 




(1) 



with the parameter S depending on the central charge c of the stable vacuum. An interesting 
question arises about the exact dependence of this exponent on c. From the exact solution 
one can figure out that 5 = 3/2 for the decay into a vacuum without massless modes, where 
c = 0. At the end of this paper we develop arguments in favor of the following exact formula 



which gives 3/2 at c = while 5 ~ (13 — c)/6 + 0(l/c) in the semiclassical limit. 

Nucleation in flat background. To begin with we recall the standard consideration 
about the decay of a metastable state in statistical mechanics, following [1], §162. Formally 
identical analysis applies to the semiclassical tunneling decay of the "false vacuum" . Here we 
restrict the discussion to the leading exponential factor, although more accurate treatment, 
which allows to estimate the pre-exponential multiplier too, can be carried out, see [3,4]. 
In [16] the two-dimensional version is specially treated in very details. 

Consider a flat infinite two-dimensional plane filled up with a metastable vacuum 4 . Its 
specific energy E\ is meant to be slightly above the energy of the "true" stable vacuum £ . 
We denote /x the (positive) difference between the two energies, i.e., E\ = So + \i. Imagine 
now a fluctuation inside this false sea, where a region T of true vacuum is formed, and let 
A(T) and p(T) be respectively the area and perimeter of T. Also, let a > be the surface 
tension at the border between true and the false vacuum. The energy of the fluctuation 
AE(T) = p(T)a — A(T)fi is a result of a competition between the gain in the bulk and the 
expense at the perimeter. Weak metastability means that the characteristic AE is much 
smaller then the temperature (in the field theoretic context AE is interpreted as the action 
of the "instanton" fluctuation and is compared with the Plank constant). In this case the 
estimate for the probability P(T) of this fluctuation is 



where T (or the Plank's constant h in the case of quantum tunneling) is the parameter 
governing the fluctuations. In what follows we always imply that the energy (action) is 
measured in the units of T (h) and omit the corresponding factors. 

By usual arguments [3,4], the specific (per unit area of the 2D space) probability of the 
false vacuum decay is controlled mainly by the configuration which extremizes Q, the decay 
instanton, a.k.a. the critical droplet. The instanton is a circle of radius 



Considerations of this paragraph go in a space of any dimensionality as well. Here we take it two- 
dimensional to match with what follows. 
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13 - c + y/(l-c)(25-c) 
12 



(2) 




(3) 



r c = a/// 



(4) 
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(figHJ. In this calculation we assume that the size of the critical droplet is large as compared 
to all microscopic scales, such as the thickness of the border of the droplet etc. This is 
certainly the case in the weak metastability limit fi —>■ (as the eq.(jl} indicates). Evaluating 
the energy of the critical droplet one finds AE C = na 2 /fi. This means that at small fi the 
decay rate is exponentially suppressed 

P~exp(-^) (5) 

in the energy gap \i. 

Nucleation in 2D Liouville gravity. Now we're going to couple our system to the 
gravitating background, where the effective action of the gravity is given by the Liouville 
action A L [15]. Choosing the coordinate system where g a b = e^&ab we have 

A L = ^-J (d a <p) 2 d 2 x (6) 

Parameter k is related to the effective number of massless modes in the corresponding vac- 
uum, called the central charge c [17]. In the semiclassical limit c — > — oo it has been found 
that k ~ — c [15]. In general the relation between k and c depends on the regularization 
prescription. Regularization with respect to a fixed background metric implies the following 
expression [18] 

k = Qb- 2 (7) 
in terms of the standard Liouville parameter b related to c through 

26- c= l + G^ + fe) 2 (8) 

As in the case of the flat background, here we are interested in the specific (per unit 
area) probability of the decay in the limit when the area is very large 5 . Correspondingly, 
our setting of the decay problem is as follows: at the beginning we have an infinite, globally 
flat space filled with the metastable phase. Semiclassically this is described by the solution 
tp = const, which we choose to be (p — 0. Imagine a local fluctuation somewhere near the 
origin, so that <p = for \x\ > r with some r. At |x| < r we have a seed of the stable 
phase, which has negative energy density — fi with respect to that of the false one. Inside 
the droplet we are looking for the extremum of the functional 

with respect to the scale function (p(x) at \x\ < r. The extremum is of course a solution to 
the positive curvature Liouville equation 

ddip = -Unfik^e^ (10) 

This limit can be achieved by fine-tuning the "bare" cosmological constant Ao to a certain critical value 
A c , so that the physical cosmological constant vanishes. That is how it is done in the matrix models of 
gravity (see e.g. example in the next section). This critical value (more precisely, — A c ) is interpreted as the 
intensive part of the specific free energy, and it is this quantity which develops the singularity /i s at /i — ► 0. 
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The general solution respecting the symmetry of the problem is a sphere 



, = 4flV , (id 

(1 + a 2 xx) 

of radius R = R c with 

* - 2k (12) 

and a scale parameter a to be determined from the continuity at the boundary of the swelling. 
It will turn instructive not to imply the relation R = R c from the very beginning considering 
/i and R as unrelated parameters, to be bound later in the extremum condition. 

We need to sew this solution with the outside region in a continuous way, i.e., at |x| = r 
we have 2Ra = 1 + a 2 r 2 . This quadric has two solutions for a 

R±VW^2 

a± = ~ 2 (13) 

Here a + corresponds to the "bubble" solution, while a_ is a "cap" one, see figO Both are 
conveniently uniformized through the parameterization 

2R = r(t + r 1 ) (14) 

In this way both branches are covered if we set ar = t and let t run from to 00, t < 1 
corresponding to the "cap" configuration and t > 1 being related to the "bubble" . It is easy 
to evaluate in this parameterization both terms of the Liouville action Q 

k f dcpdifd 2 x = \ (log(l + t 2 ) - (15) 



24tt J r r 6 V 1 + 1 2 

H I e Lp d 2 x = ^\iR 2 



x\<r 



1+t 2 



Adding the surface tension term 2nar we find the total energy of the swelling to be varied 
E= k - (log (1 + t 2 ) - ^ - x 2 {\ + t 2 ) + 2s*) (16) 
were we have used the notation 



•W£ (17) 



instead of the surface tension. Also we introduced another parameter 

2x = r/R c . 



6 



The relevant extremum is achieved at 



t = s and x - 
value 

E r 



8/(1 
k 



The extremal 
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log (1 + s 2 ) 



(19) 



It can be verified directly that in both cases 
the chosen extremum corresponds to a sad- 
dle point in the parameters (R,r), in the 
case s <C 1 the unstable direction being 
mostly along the "perimeter" direction r, 
while in the bubble case s > 1 it is generally 
along the "inflation" direction R. Finally, 
we find for the decay rate 



1 + 



67RX' 

kfj, 



-k/6 



(20) 




Figure 3: 
sphere g - 
triangles. 



An example of triangulation of a 
0. Ising spins are attached to the 



At the extreme limit k — > oo we are back 
to the classic droplet result (jSJ. On the 
contrary, the weak metastability limit \i — ► 
is characterized by a power-like suppression P ~ (k^i/ a 2 ) k ^ . 

Gravitational Ising vacuum energy. We recall in this section very briefly the dy- 
namical lattice Ising model (DLIM), as it has been proposed in [14], and present without 
derivation the form of its exact solution in the scaling limit near the ferromagnetic phase 
transition. The formulation is standard for the two dimensional random lattice models and 
starts with an ensemble j of planar graphs (in general of genus g) and of size N. It is 

considered well established that the scaling singularity, which is supposed to correspond to a 
continuous limit, doesn't depend on the details of this ensemble (neither on most microscopic 
details of the spin model located on the graph). Here for definiteness we will imply a very 



particular realization of 

tions, like an example in fig 
number of triangles in G^ 9: 




as the graphs made exclusively of triangles, the triangula- 
In this case the size N can be given precise meaning, the 



N 



In DLIM we associate a spin variable Oi = ±1 with each triangle i = 1, 2, . . . , N. As in the 



usual Ising, only "nearest neighbor" triangles (ij) 
to spin-spin interaction 



(i.e., having a common edge) contribute 



w[W] 



EK<Ji<7j 
(ij) 3 



(21) 



where K is the "thermal" parameter (called traditionally the exchange integral) and H is the 
external magnetic field. In this study we consider the large iV limit characteristics, namely 
the specific free energy (per triangle) of DLIM. This energy is an intensive characteristic and 
as such doesn't depend of global properties, like g. For our purposes it is most convenient to 
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start with the spherical partition function (g = 0), where the result of [14] applies directly. 
Introduce the microcanonic spherical partition function 

Z N (K,H) = J2f G ^J2 {rT} ex V ( <- n ti (T W (22) 

Here, in addition to the standard Ising sum over the spin configurations {a} the sum over 
| is also implied, the last being the lattice version of the quantum gravity path integral 
over the metrics. Considering the leading logarithmic behavior at large N 

log Z N (K, H) ~ -N£(K, H) (23) 

we readily see that the intensive energy S(K, H) is determined by the position of the right- 
most singularity in the grand partition function 

Z(K, H, X) = J2 N e~ XN Z N (K, H) (24) 

in the "chemical potential" X. 

Let (X, K, H) = (X c , K c , 0) be the position of the double singularity in the grand parti- 
tion function, caused by the combined divergence in size N in (|2^jl and of (|2^|l due to the long 
correlations of spins at the ferromagnetic phase transition. The singular part Z sing (x, t, h) of 
the partition function depends in a scaling way on the deviations x ~ X — X c , t ~ K c — K 
and h ~ H (precise definitions depend on the choice of scale). Boulatov and Kazakov [14] 
showed that the function 

u = d 2 Z sing /dx 2 (25) 
solves the following simple algebraic equation 

x=u3+ h 2+ wU (26) 

As usual the continuous theory deals with the singular part of the intensive energy, which, 
up to sign coincides with x c , the singularity of equation (|2l)|) . It is easy to see [13] that the 
relevant singularity is located at 

x c = |/ (5/ 2 + 9/ + 3) (27) 

where / is an appropriate solution to the following fifth order algebraic equation 

2h 2 

f(f + l) 4 = W (28) 

As it is argued in ref. [13], the low temperature phase t < of DLIM corresponds to the 
solution 

A r(W4-i) / 2V^ y 



and leads to the following singularity at h — > 



t x c 



1 oo 

+ 3 £ 



n 



l)r(5n/4 - 3) / 2 1 / 4 / l V2 



n=2 

1 VQh 

2 ~ (-t) 5 / 2 + 3 3 / 4 (-t) 15 /4 



n!r(n/4) 

23/4^3/2 



+ 



3 l/4(_ t )5/4 
/i 2 



(30) 



4(-t)< 



+ 



The free energy exhibits the square root branch cut at h — with the exponent 3/2. 

Discussion and Speculations. There are serious reasons to believe that in two di- 
mensions the Liouville gravity affects the standard nucleation mechanism of the metastable 
decay in a significant way. The standard exponential suppression of the decay rate is substi- 
tuted by a weaker power-like one, changing qualitatively the picture. We were only able to 
obtain more precise information either in the semiclassical limit, controlled by the instanton, 
or in the exactly solvable dynamical lattice Ising model. One could expect in principle a 
one-parameter family of universality classes, depending on the central charge in the final 
state of the decay. A better understanding of this family is desired. In this relation dif- 
ferent steps would be valuable. One thing, which is not attempted in the present study, 
is the one-loop (and higher) corrections to the critical swelling instanton. Such calculation 
would provide a constant term in the semiclassical expansion 5 ~ — c/6 + const of the ex- 
ponent in (JTJ), and thus support (or disprove) the conjectured exact form (J2J). It is also 
useful to look for more solvable examples among the dynamical lattice statistical systems 
exhibiting metastability and, preferably, having non-trivial massless modes in the stable vac- 
uum. And of course, metastability opens a large playground for Monte-Carlo simulations, 
essentially dynamical processes. Many appropriate models are known (e.g., the "gener- 
alized random lattice Ising model of ref. [19]), which certainly can be put to metastable 
state, while the central charge in the stable vacuum remains completely under our control. 

The Eq.((T} applies to the decay 
probability in the situation when the 
metastable phase is accommodated by 
infinite globally flat geometry, or more 
generally when all the length scales as- 
sociated with that geometry are large 
as compared with In particular, 

the physical cosmological constant of the 
metastable phase must be much smaller 
then \i. This settlement is perfectly ap- 
propriate for the situation of the previous 
sections, where we were interested in the 
intensive part of the specific free energy. 
But the semiclassical calculations of the third section certainly can be generalized to incor- 
porate the effect of the cosmological constant in an interesting way. 




Figure 4: Critical swelling in the regime of ex- 
tremely weak metastability fi — > 0. 
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To conclude, we speculate a conjecture for the exact relation (J2J) between the central 
charge and the exponent in the decay law. It is easy to see, that in the limit of extremely 
weak metastability /i — > the "bubble" configuration of fig. Efb) tends to an almost complete 
sphere, connected to the metastable plane through a small "throat" . The last can be seen as 
a pointlike puncture (figEJ). We recall in this context that the partition function of a compact 
sphere behaves as fi b +1 where b is from (jHJ). This is an exact result of quantum Liouville 
field theory, which can be found e.g., in [20]. A puncture corresponds to a differentiation 
with respect to [i and the action of the limiting swelling is evaluated as /x 1 ^ , i.e., the power 
law conjectured in (J2J). In this consideration we assume (with no really good reasons) that 
the essentially semiclassical arguments about the critical swelling, which we used to reduce 
the problem to a punctured sphere, apply qualitatively also at general quantum c. 

Let us mention also that the exponent equal to 1/b 2 + 1 (which differs by 1 from (J2J)) 
was proposed recently by A.Polyakov [21] in apparently different context, in relation with 
possible unstability of the de Sitter space. Obviously, understanding precise relation to the 
decay problem considered here would be of much interest. 
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